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ABSTRACT 

Multi-zone models of Type I X-ray bursts are presented that use an adaptive nuclear reaction network 
of unprecedented size, up to 1300 isotopes, for energy generation and include the most recent mea- 
surements and estimates of critical nuclear physics. Convection and radiation transport are included 
in calculations that carefully follow the changing composition in the accreted layer, both during the 
bursts themselves and in their ashes. Sequences of bursts, up to 15 in one case, are followed for two 
choices of accretion rate and metallicity, up to the point where quasi-steady state is achieved. For 
M = 1.75xl0~ 9 M© yr _1 (and M = 3.5xl0 -10 M yr _1 , for low metallicity), combined hydrogen- 
helium flashes occur. These bursts have light curves with slow rise times (seconds) and long tails. 
The rise times, shapes, and tails of these light curves are sensitive to the efficiency of nuclear burning 
at various waiting points along the rp-process path and these sensitivities arc explored. Each displays 
"compositional inertia" in that its properties are sensitive to the fact that accretion occurs onto the 
ashes of previous bursts which contain left-over hydrogen, helium and CNO nuclei. This acts to re- 
duce the sensitivity of burst properties to metallicity. Only the first anomalous burst in one model 
produces nuclei as heavy as A = 100. For the present choice of nuclear physics and accretion rates, 
other bursts and models make chiefly nuclei with A 64. The amount of carbon remaining after 
hydrogen-helium bursts is typically < 1% by mass, and decreases further as the ashes are periodically 
heated by subsequent bursts. For M = 3.5 xlO -10 M© yr _1 and solar metallicity, bursts are ignited in 
a hydrogen-free helium layer. At the base of this layer, up to 90% of the helium has already burned to 
carbon prior to the unstable ignition of the helium shell. These helium-ignited bursts have a) briefer, 
brighter light curves with shorter tails; b) very rapid rise times (< 0.1s); and c) ashes lighter than 
the iron group. 

Subject headings: neutron stars, X-ray bursts, nucleosynthesis 



1. INTRODUCTION 

Since the 1970's, when the nuclea r instability of accre- 
tion onto neutron stars was noted ({Hansen fc van Horn! 
Il975|) . and the connection with short, transient 
X-ray flashes (Type I X-ray bursts) pointed out 
ijWooslev fc Taaml |1976|L there have been numerous 
studies of the rmonuclear flashes on neutron star s . For 
reviews, see IBildstenl 17553 iLewin et al.1 EH EH 
iStrohmaver fc Bildster]l2003l Studies in the mid-1980's 
elucidated the relevant nuclear phys ics, known as the 
rp-process (|Wallace fc Wooslevlll981|h and showed that 
the most critical parameter determining burst proper- 
ties was the accretion rate, with s everal different regimes 
of bursting behavior expected ijFuiimoto et al.1 119811 
iFushiki fc Lamb|ll987|) . In particular, combined helium 
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and hydrogen burning powers flashes with the lowest 
critical masses for accretion rates between 4-11 xlO -10 
and 2xl0 -8 M Q yr _1 (depending upon the metallicity 
of the accreted matter) , whereas pure helium flashes oc- 
cur beneath a stable hydrogen shell for rates between 
2xl0- 10 and 4-11 xlO~ 10 M Q yr _1 . Weakly flashing hy- 
drogen shells occur for rates less than 2xl0 -10 , M Q yr _1 . 
The generic properties predicted for these flashes — in- 
tervals, durations, energies, brightness, etc. — were in 
good agreement with observations, though the trends in 
burst properties with changing accretio n rate were of- 
ten at odds with the sim ple theory (e.e.. lBildsterJl2(T0fi: 
Ivan Paradiis et al.lll98l . 

Previous theoretical studies can generally be criti- 
cized, however, either for oversimplification of the nu- 
clear physics, especially the use of small approximation 
networks for the e nergy generation, e.g.. fAyasli fc .Tossl 
Il982l iTaanl UMi IWooslev fc Weaver! 11984 or for in- 
adequate attention to the stellar model, especially con- 
vection, i n single-zoned studies that use large rea ction 
networks ijBrown et al.ll2002HSchatz et alJl2001albTl . In 
this paper, we combine both: large networks and current 
nuclear data with leading edge (albeit one-dimensional) 
stellar models that include convection and multi-zone ra- 
diation transport. 

Any such attempt to simulate Type I X-ray bursts real- 
istically meets with four challenges. First is the physics 
of the accretion process and the geometry of the run- 
away. Is the accretion uniformly distributed over the 
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surface of the neutron star prior to runaway and does 
the burst commence almost simultaneously over that sur- 
face? Observations of nearly-coherent oscillations dur- 
ing Type I bursts in the last six years have brought 
these questions into focus (see iStrohmaver fc Bildstenl 
120031 for a review). The oscillations, which are inter- 
preted as rotational modulation of brightness asymme- 
tries, indicate rapid rotation, and perhaps non-uniform 
ignition. If so, rotation may be key to balancing the 
strong trans verse pressure gradients that arise near a lo - 
cal hot spot <|Spitkovskv et al.ll2002tlZingale et al.ll2003|) . 
We do not address these issues in this paper; our one- 
dimensional calculations address only the local physics 
per unit area. While comparison to observables like the 
light curve could be affected by the geometry, the nuclear 
physics of the runaway is not. 

Second, the nuclear physics is complex with no single 
or even several reactions governing the energy genera- 
tion rate. Recent years have brought significant advances 
in our understanding of the major flows in the rp- and 
ap-proccsses a nd the properties of the nuclei that goy - 
ern them fe.g.. lBrown et al.l2002USchatz et al.l2001albfl . 
These advances should be included in any modern study. 

A particularly exciting development has been the dis- 
covery of very energetic, long duration Type I bur sts, 
known as "superbursts" (see IKuulkers et alJ 120021 for 
a review of supcrburst properties). These flashes 
have 1000 times the energy and duration of normal 
Type I bursts, and have been in terpreted as unstable 
ignition of a thick carbon layer ijCumming fc Bildstenl 
| 200lHStrohmaver fc Brownl2002|) . originally proposed by 
IWooslev fe Taaml (|1976j T as a gamma-ray burst model. 
Calculating the amount of carbon produced during un- 
stable hydrogen/helium burning, and how much carbon 
survives to great depth, requires an accurate multi-zone 
calculation with a large nuclear network. 

The third challenge is following not just o n e, bu t 
many bursts. As was pointed out by iTaaml l)1980j) . 
the thermal and compositional "inertia" of the neu- 
tron star envelope can have important implications for 
the properties of subsequent bursts. It is only by 
computing a succession of X-ray bursts that the heat- 
ing associated with the previous thermonuclear flashes 
and compositional ch ange in the accreted layer can 
be taken into account |Avasli fc Josslll982t iTaam et alJ 
119931 IWooslev fc Weaverlll984|) . Here we not only carry 
out fine-zoned studies of individual bursts, but follow 
each case until a steady state is achieved, including in 
one study, 15 consecutive bursts. 

The fourth challenge is making contact with the rich 
archive of photometry and spectra for observed X-ray 
bursts. The models presented here are limited to single- 
temperature black bodies calculated using flux-limited 
radiative diffusion, although more detailed studies of the 
spectrum and color temperature could use our results as 
input to a more superior treatment of radiation trans- 
port. An immediate application of our results is direct 
comparisons with observed light curves to study, for ex- 
ample, the physics underlying burst rise and decay times. 
We do not attempt such a comparison in this initial pa- 
per, but leave this for future work. 

An outline of the paper is as follows. In § El 
we describe modifications to the onc-dimcnsional im- 
plicit hydrodynamics code - Kepler - used for this 



study. This code has been used before to study X-ray 
bursts re.g..lFushiki et al.llT99l lTaa,m et al .llT9M ITflM 
IWooslev fc Weaveill984|) and we need only clarify recent 
modifications to the nuclear reaction library and the im- 
plementation of energy generation from an extended net- 
work. A novel approach to evolving the abundance vec- 
tor in each zone, called an "adaptive reaction network" , 
is employed in which the number of isotopes employed in 
each cycle of each stellar model may vary according to 
the active flows and significant abundances present. At 
maximum, we employed over 1300 isotopes in the reac- 
tion network for over 1000 zones at a time. The energy 
generation obtained from this network is used in calcu- 
lating the stellar model. 

In subsequent sections we describe the results for four 
models that crudely sample the accretion rates and 
compositions relevant for X-ray bursts. Two accretion 
rates, 3.5 and 17.5 xlO -10 M Q yr _1 and two metallici- 
ties, Z = 0.001 and 0.02 are examined. As expected, 
in three of these combined hydrogen helium runaways 
were observed. In the fourth, Z = 0.02 and M = 
3.5xl0 -10 M Q yr _1 , hydrogen burns away before helium 
ignites in a thick shell that has already burned, at its 
base, mostly to carbon. 

In the results and conclusions sections we describe 
some of the novel aspects of our results, the sensitiv- 
ity to nuclear rates employed, the light curves and in- 
tervals for the bursts, and the composition of the ashes 
after many bursts have occurred, addressing both the 
expected composition of the crust, and whether enough 
carbon is produced to power a superburst. 

2. THE KEPLER CODE AND THE NUCLEAR DATA 
EMPLOYED 

Kepler is a one-dimensional implicit hydrodynamics 
code adapted to the study of st ellar evolution and ex plo- 
sive astrophysical phenomena ((Weaver et al.lll978l) . It 
has been used for twenty years to simulate X-ray bursts 
on neutron stars. The equation of state allows for a gen- 
eral mixture of radiation, ions, and electrons of arbitrary 
degeneracy and relativity Electron-positron pair forma- 
tion is also accurately included. Convective mixing is 
modeled using mixing length theory in a time-dependent 
manner. That is, the composition is diffusively mixed 
with a diffusion coefficient calculated from the convec- 
tive velocity. The convective criterion is Ledoux but 
with a substantial semi-convective diffusion coefficient, 
about 10% of the radiative diffusion coefficient, in re- 
gions that are stable by the Ledoux criterion and un- 
stable by Schwarzschild. A relatively small amount of 
convective overshoot is included by flagging convective 
boundary zones as semi-convective. The opacity, neu- 
trino loss es, and other aspects of the code have been dis- 
cussed in lWooslev et all l|2002j) . In particular, the recent 
OPAL and Los Alamos opacity tables are used wherever 
the helium mass fraction exceeds 10 _J and the tempera- 
ture is less than 10 9 K. 

A ccretion is handled as in IWooslev fc Weaver! l(1984f) 
and ITaam et alJ l|1996th In the Lagrangian code, the 
outer boundary pressure is increased to simulate the 
weight of the accreted matter at the given gravitational 
potential. This continues until the pressure equivalent to 
the mass of a new zone has accumulated. At that point a 
new zone is added to the grid with conditions like those in 
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the previous outer zone. The surface boundary pressure 
is zeroed and the process begun anew. All runs accrete 
zones of 2xl0 19 g and all accreted zones are retained un- 
til the end of the run, i.e., there was no re-zoning of the 
accreted layers. Some models accumulated more than 
1000 zones this way. Typically, the zones at the base 
of the accreted layer become significantly thinner than 
1 cm! 

During the accretion process, nuclear energy genera- 
tion and composition change are calculated at every time 
step. Convection is always "on" in the sense that, if the 
Ledoux criterion for instability is satisfied, the code re- 
sponds. In this way we found that, even though the 
strongest convection accompanied the onset of a burst, 
there were interesting episodes of both convection and 
thcrmohalinc mixing in between the bursts. 

2.1. Adaptive network 

The major improvement in the present study was a 
more detailed treatment of the nuclear physics and en- 
ergy generation. Unlike previous multi-zone calculations 
that carried only a very sparse "approximation network" , 
we computed energy generation and composition changes 
with a much larger network previously used only to study 
nucleosynthesis. 

Reaction rates were extracted from a library of nuclear 
data (§ 12.3(1 carried in the calculation for isotopes rang- 
ing from hydrogen to polonium (Z = 1 to 84). When 
the abundance of an isotope exceeded a threshold mass 
fraction of 10 -17 , all neighboring isotopes that could be 
accessed by (n/y), (p,7), (a/y), (a,n), (a,p), or (p,n) 
reactions and their inverse, as well as (2p,7) were added. 
For (n,7) and (7,n) links the limit is somewhat smaller, 
10~ 18 . Additionally, all possible decays, (e~,v c ), (e + v c ) 1 
{e~v e ) 1 and (a), were included for all isotopes in the net- 
work independent of their abundance. Conversely, when 
the mass fraction of an isotope dropped below 10~ 23 , 
and its presence is not necessary to satisfy any of the 
above conditions for adding it, it was removed from the 
network. These criteria were applied in all zones during 
each stage of the calculation and only one network, the 
sum of all these local conditions, was used at any point in 
time throughout the star. This is necessary since zones 
may mix at unpredictable times and locations during the 
evolution. 

2.2. Binding energies 

Nuclear flow through the rp-process waiting point nu- 
clei is exponentially sensitive to the proton separation 
energies (S p ) of a few nuclei near the proton-drip line. 
Expe rimental values fo r the b inding energies were taken 
from lAudi fc Wanstral <)1995t AW95) where available 
(Fig. For many of the relevant nuclei however, no 
experimental information is available and one must rely 
on theoretical mass predictions. For nuclei with Z > N 
in the mass range A = 41-75, we used the compilation of 
iBrown et alJ l|2002(L in which the experimentally deter- 
mined masses of N > Z nuclei, together with a Skyrme 
Hartree-Fock calculation of Coulomb displacement ener- 
gies, are used to estimate masses for the Z > N mirrors. 
These mass predictions have an uncertainty of about 100 
keV for the mass difference to the mirror nucleus. This 
results in an average uncertainty of the proton separa- 
tion energies of about 140 keV if the mirror masses are 



sufficiently well known. As Table ^ shows, for many of 
the relevant rp-process nuclei this is not the case and un- 
certainties are substantially larger, though still within a 
few hundred keV. 

For another large s et of nuclei, we used the unpub- 
lished calculations of <)Brownll2002fl (Fig. 0J. These are 
obtai ned using the same displacement energy method 
as in IBrown et alJ l|2002|) . but the theoretical error in 
the displacement energy is larger than 100 keV (though 
still within a few 100 keV) because of the need to ap- 
ply a spherical basis to a mass region where some nuclei 
are strongly deformed. Furthermore, for many of these 
heavier nuclei, the masses of the N > Z mirrors are 
not known experimentally. We then used the theoretical 
extrapolations of AW95 as a basis for the displacement 
energy method. The errors in the AW95 extrapolations 
for A > 76 introduced an additional error of order 500 
keV for nuclei near the Z = N line. When no information 
is available either from AW95 or Brown, we employed the 
theoretical estimates of iMoller et al.l (|1995l) . 

Table ^ shows the proton separation energies of sev- 
eral key nuclei important in determining flow past the 
waiting point nuclei with long half-lives against positron 
decay - 6 4 Ge 68 Se, 72 Kr, and 76 Sr. A s discussed by 
lOrmandl (|1997j) and IBrown et all J2002) , uncertainties 
in the masses of 64 Ge and 68 Se dominate uncertainties 
in the proton separation energies of 65 As and 69 Br, key 
nuclei in the rp-process. 

2.3. Reaction rates 

Nuclear reaction rates were taken from experiment, 
shell-model calculations, and statistical model (Hauscr- 
Feshbach) calculations. An experimentally determined 
rate was adopted whenever possible 7 . Largely, though, 
experimental information was unavailable for the proton- 
rich nuclei important in X-ray bursts. For nuclei in the 
mass range A — 44-63, experimentally-undetermined 
(p,7) ra tes were tak e n fro m the shell model calcula- 
tions of iFisker et alJ (|£001). Those authors also pro- 
vide a critical discussion of uncertainties inherent in the 
different methods of calculating rates near the proton 
drip line. All other reaction rates were calculated using 
the Hauser-Feshbach code Non -Smoker as described by 
iRauscher fc Thielemann! <|2000l) . 

2.4. Weak rates 
Our calculations include electron capture 

er + A Z^ A (Z-1) + ^ C , (1) 

nuclear e + decay (positron emission) 

A Z^ A (Z-l) + e+ + ^ e , (2) 

and the neutrino energy loss rates associated with the 
above processes. For low temperature and density (Tg < 
0.5, pY e < 10 5 gcm -3 ), we adopt the ground state rates. 
This is appropriate for the proton rich nuclei of inter- 
est during rp-process burning. These nuclei have decays 
characterized by large Q- values, for which positron emis- 
sion dominates over electron capture. We take ground 
state weak decay rates from experiment where available 

7 See http://www-pat.llnl.gov/Research/RRSN/ for sources of 
experimentally determined rates. 
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and from the comp ilatio n of IMoller et all l)1997|) other- 
wise. IMoller et al.1 <)1997l) calculated weak decay assum- 
ing only the presence of Gamow- Teller (GT) transitions. 
As is discussed below, Fermi transitions dominate the 
lifetimes of many near-proton-drip line nuclei. For such 
nuclei the lifetime is typically a factor of 4-6 shorter than 
that estimated from a consideration of GT transitions 
alone. This issue and others relating to the weak physics 
of TV ~ Z nuclei will be addressed in a future study. 

For nuclei with A < 60 and for higher temperature 
and density (T 9 > 0.5, pY e > 10 5 gcm~ 3 ), we in- 
clude the influence on the weak rates of thermal pop- 
ulation of excited states as well as continuum electron 
ca pture. These rates are taken from the compilation 
of lLanganke fe Martinez-Pinedol J200C , LMP), and from 
the compilatio n of iFuller et alJ 111982 , FFN) for nuclei 
not studied bv lLanganke fc Martinez-Pinedol With the 
exception of a slightly different interpolation of neutrino 
energy loss rates, interpolati on of the we a k rate s is done 
following the prescription of IFuller et al.l l|1985|) . 

Thermal effects on weak rates for near-proton-drip line 
nuclei in the mass range A ~ 60-100 follow simple sys- 
tematics. Consideration of these systematics can be used 
to determine if our use of ground state lifetimes for these 
nuclei is reasonable. Fig. [2] illustrates important nu- 
clear and weak flows near the waiting point nucleus Kr. 
These flows arc typical of those near the other waiting 
point nuclei in this mass region. 

As can be seen from Fig. El important weak decay par- 
ents fall into four categories. Nuclei with Z = TV + 1 
undergo e + -decay to their isospin mirrors. For these nu- 
clei, the ground state of the (e + ) daughter is the isobaric 
analog state (IAS) of the parent ground state, the first 
excited state of the daughter is the IAS of the first excited 
state in the parent, and so on for all of the levels. Because 
the nuclear Hamiltonian commutes with the isospin rais- 
ing and lowering operators, the Q-values and rates for all 
of these (parent level) — > (IAS in daughter) Fermi transi- 
tions are approximately the same. GT transitions are 
not expected to dominate the decay rates because the 
GT strength is typically spread out over a wide range 
in daughter excitation energy. In addition, electron cap- 
ture cannot compete with these large Q-value (3 + decays 
at the small temperatures and electron Fermi energies 
achieved in X-ray bursts. Consequently, the decay rates 
of Z = TV + 1 nuclei are essentially temperature and 
density independent (for the range of conditions found 
in X-ray bursts). For nuclei with Z = TV + 2 = even, 
the situation is similar. Each parent state decays via a 
Fermi transition to an IAS in the daughter and the de- 
cay rates for these nuclei are well approximated by the 
ground state decay rate. 

Nuclei with TV = Z = odd also decay principally via 
Fermi transitions. However, while for these nuclei it is 
true that every low-lying daughter state has a low-lying 
IAS in the parent, the converse is not true. Parent states 
in these N = Z = odd nuclei that do not decay via a 
Fermi transition typically decay about an order of mag- 
nitude slower than do states with a low lying IAS in the 
even-even daughter. Consequently, e + -rates for these nu- 
clei are approximately proportional to Gdaughtcr/G pare nt 
l)Pruet fc Fuller! 120021) . Here G is the partition function 
and the ratio is a rapidly decreasing function of temper- 
ature because of the high level density of odd-odd nuclei 



(relative to even-even nuclei). 

The last, and most important, class of Z ~ TV nuclei 
are those with Z = TV = even. Of these, only two have 
partition functions at Tg = 1.5 ( 76 Sr with G ss 1.7 and 
80 Zr with G « 1.6) large enough for the de cay of excited 
states to determine the lifetime (see also iSchatz et alJ 
l)1998(l for a discussion of 2 + lifetimes in X-ray bursts and 
estimates of 2 + excitation energies for Z > 82 nuclei). 
For these two nuclei the thermal decay rate may be larger 
by a factor of ~ 2 — 3 than the ground state decay rate, 
though a more reasonable estimate is probably in the 
30 % range. As with the other important decays, electron 
capture is negligible for TV = Z = even nuclei during rp- 
process burning. 

Detailed calculations of thermal effects on the weak de- 
cays of A > 60 nuclei should eventually be incorporated 
into X-ray burst studies. However, the above consid- 
erations suggest that the adoption of laboratory ground 
state rates is not too bad. One shortcoming of our proce- 
dure concerns TV = Z = odd nuclei. We typically under- 
estimate the lifetime of these nuclei in the X-ray burst 
environment by about a factor of 5. For A < 80 this 
error is not expected to be important because so little 
of the nuclear flow (~ 1 % or less) proceeds through the 
e + -decay of these nuclei. With increasing atomic mass 
number the proton drip line tends closer to TV = Z . Er- 
rors in the lifetimes of TV = Z = odd nuclei may have 
a larger influence on nuclear flow for A > 80. However, 
a qualitative change in the flow is unlikely because the 
typical lifetimes of important A > 80 nuclei are so short 
(~ 1 s) compared to lifetimes of long-lived waiting point 
nuclei (« 64 s for 64 Gc). 

2.5. Neutrino Losses from Radioactive Decays 

Neutrino losses during weak decays are an important 
part of the total energy budget, since neutrinos typically 
carry away a sizeable fraction (half or more) of the en- 
ergy available in a decay. When a weak decay rate is 
estimated from the compilation of LMP or FFN, we also 
adopt the LMP or FFN value for the neutrino energy 
loss rate associated with that decay. For a few light nu- 
clei ( 7 Be, 13 N, 14 ' 15 0, 17 ' 18 F), neutrino energy loss rates 
are calculated using experimentally determined ground 
state weak strength distributions. Other neutrino en- 
ergy loss rates were estimated using a code provided to 
us by Petr Vogel. This code uses an empirical form of 
the GT + strength distribution to estimate the average 
energy of neutrinos emitted in a decay. The neutrino en- 
ergy loss rate is then taken to be the product of the weak 
decay rate and the average neutrino energy Though em- 
pirical strength distributions cannot reliably reproduce 
ground state decay rates, they can do a fair job of es- 
timating average neutrino energies because phase space 
so heavily favors those transitions with the most ener- 
getic outgoing neutrinos. For decays characterized by 
large Q-values, Vogel's estimates of average neutrino en- 
ergies typically agree with more detailed estimates (FFN, 
LMP) to within 10% or 15%. 

3. INITIAL MODELS AND SETUP 

Four different combinations of accretion rate and 
metallicity were examined (Table GJ- The first 
model, zM, employs conditions very similar to 
ISchatz et al.l <)2001a|k an accretion rate of O.lMEdd = 
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8.8xl0 3 gcm -2 s -1 = 1.75xlO -9 M yr -1 for a 10km 
radius neutron star. For comparison, we also calculated 
a similar model, zm, with an accretion rate one-fifth as 
large. The composition in both cases was taken to be 
75.9% hydrogen, 24% helium, and 0.1% 14 N by mass, 
correspondin g approximately to the Z = 10 -3 used by 
iSchatz et al.L Hereafter this is referred to as the 5 % so- 
lar composition, since the total mass fractions of CNO 
in the sun are about 20 times larger. For comparison, 
we also calculated similar models ZM and Zm with solar 
metallicity (70.48% 1 H, 27.52% 4 He, 2 % 14 N) and high 
an d low accretion rates respec tively. These are similar to 
the lAnders fc Orevessel (|1989ft solar abundances, but the 
rearrangement of the CNO isotopes into nitrogen that 
naturally occurs in the CNO cycle was skipped since this 
will occur rapidly at the high temperatures of the ac- 
cretion process and the energy added is inconsequential 
compared with that of accretion. 

To save time and improve energy conservation, only 
the outer 2xl0 25 g of neutron star crust is carried in the 
calculation. This is orders of magnitude larger than the 
mass of all X-ray bursts combined, so the layer essen- 
tially acts as a large repository of thermal inertia. Its 
composition is taken to be iron and no nuclear reactions 
are followed in this layer. We take the luminosity at 
the base of the substrate to be 3.2 xlO 33 ergs -1 for Mod- 
els zm and Zm and 1.6xl0 34 ergs -1 for Models zM and 
Z M, i.e., the ac cretion rate times 0.15 MeV/nucleon as 
in ISchatz et all Before accretion is switched on, we al- 
low the substrate to relax to thermal steady state with 
the power input at its base balanced by the luminosity 
flowing into the accreted zones. 

A 1.4M neutron star with radius R = 10km is 
adopted giving a surface gravity g = GM/R 2 = 
1.9 xlO 14 cm -2 s -1 . We do not include the effects of gen- 
eral relativity in our simulations, we address the rela- 
tivistic corrections that must be applied to our results in 

§1331 

4. RESULTS 

4.1. Model zM 

The model considered in greatest detail and followed 
though the largest number of flashes was Model zM, 
w hose parameters duplic ated those of previous studies 
bv lSchatz et alJ l)2001albfl . The accretion rate and metal- 
licity (Table imply that hydrogen will survive to the 
depth where helium ignites, so that this should be a rep- 
resentative case of a combined hydrogen-helium runaway 
(§0- Previous studies have also suggested that the nu- 
clear flows should lead to the creation of quite heavy 
nuclei, with the rp-process terminating in the SbSnTc 
cycle, making this an interesting test case for the large 
adaptive network. 

4.1.1. The first burst in Model zM 

Following accretion for 41,640 s, during which a layer 
of 4.66 xlO 21 g accumulated, the temperature at the base 
of the accreted material began to rise rapidly compared 
with the accretion time scale. The hydrogen and he- 
lium mass fractions at the base of the accreted layer at 
that sample time (Fig. were 0.693 (down from an ini- 
tial 0.759 due to stable hydrogen burning) and 0.283. 
The temperature there was 2.67xl0 8 K, and the den- 
sity, 8.87xl0 5 gcm -3 . Most of the energy generation at 



this point, 3.5xl0 35 ergs -1 , was still coming from the (3- 
limited CNO cycle, though helium burning had increased 
the abundance of catalytic nuclei to mass fractions 1.19 % 
of 14 O and 1.29% 15 O, respectively. 

Fifty-three seconds later, at 41,693 s, the temperature 
at the base of the hydrogen had increased to 3.27 xlO 8 K 
and the maximum was starting to move outwards. The 
energy generation had risen by a factor of 5. Six seconds 
later, at 41,699 s, energy began to be transported by 
convection (Fig.^J when the maximum temperature was 
3.97xl0 8 K. This convection began 5xl0 20 g, or about 
half a meter above the base of the hydrogen layer. At this 
point, the mass fractions at the base of the accreted layer 
were l K: 0.683, 4 He: 0.256, 14 0: 0.037, and 15 0: 0.016. 
Over the next 7 seconds the convective region grew to 
encompass the entire accreted layer, with the exception 
of the outer 2 or 3 zones and a mass of ~ 5xl0 19 g. As 
the convection neared the surface the observable tran- 
sient commenced (Fig.^J. At this time 41, 705 s (Fig.EJ), 
the composition still consisted almost entirely of hydro- 
gen and helium though appreciable amounts of heavy 
elements were beginning to be synthesized. The nu- 
clear power being generated was 4.5xl0 38 ergs" 1 , but 
only a small fraction of this had appeared at the surface 
(L = 1.5xl0 35 ergs" 1 ). Convection ceased during the 
next second. 

From that point on the luminosity rose slowly, by dif- 
fusion, to nearly the Eddington value, 2xl0 38 ergs -1 . 
Qualitatively, it seems that the energy is transported 
by convection until an adiabatic temperature gradient is 
established to the surface. Accomplishing this requires 
raising the temperature which leads to expansion. The 
necessary PdV work against the enormous gravitational 
potential of the neutron star uses up most of the nu- 
clear energy release in the early stages of the run away. 
Once this has been accomplished, convection shuts off 
and is unimportant during the remainder of the burst. 
The light curve, for this model, is powered by diffusion in 
near steady state with the nuclear power. As the burning 
region cools gravitational potential energy is converted 
back into heat, but this is a small fraction of nuclear 
energy until the burst is essentially over. 

At 41,715 s, for example, about 5-10 s into the burst, 
the luminosity was 1.5 xlO 38 ergs -1 and all convection 
had ceased. Nuclear energy generation, 1.6xl0 38 ergs -1 , 
was in near steady state with the luminosity. An appre- 
ciable gradient in hydrogen was beginning to develop. At 
the bottom of the layer, the most abundant mass frac- 
tions were X H (0.279), 4 He (0.0 77), 60 Zn (0.460), and 
64 Ge (0.116), but halfway out (in mass) the most abun- 
dant mass fractions were : H (0.593), 4 He (0.175), 56 Ni 
(0.030), and 60 Zn (0.156). The temperature and den- 
sity at the bottom of the layer were 1.51xl0 9 K and 
4.37xl0 5 gem -3 , near the maximum developed during 
the burst. 

Energy generation from the rp-process and transport 
by radiative diffusion continued to power a brilliant dis- 
play with a long tail lasting several hundred seconds 
(Fig. [SJ) (n.b., the luminosity from the burst must be 
added to the accretion luminosity, 2xl0 37 crgs -1 here, 
in all plots of the light curve). The composition 290s 
after the burst began (at 41,992 s) is shown in Fig. 
At this point energy is continuing to be generated by 
nuclear reactions happening in unburncd hydrogen well 



6 



Woosley et al. 



above the base of the accreted layer — where hydrogen 
has burned away. Energy generation from the decay of 
radioactivities produced at the bottom of the layer is 
negligible. The temperature and density of the layer at 
this time are shown in Fig. |H1 

4.1.2. The second burst in Model zM 

Following the first burst, accretion continued at 
1.75xl0 -9 M Q yr _1 until a second critical mass had ac- 
cumulated (Table HJ. However, this time accretion oc- 
curred onto the ashes of the previous burst, which con- 
tained unburned hydrogen and helium, rather than onto 
an inert substrate. This distinction greatly altered the 
conditions for, and nature of, the second and subsequent 
bursts. 

Fig.EJshows the major abundances at 54, 592 s after ac- 
cretion began, about 30 s before the second runaway. Un- 
burned hydrogen and helium are abundant in the outer 
ashes of burst 1 - as they must be in any burst that 
has not remained fully convective throughout its burn- 
ing. Not only are hydrogen and helium abundant, but so 
are 14 and 15 O from helium burning during the inter- 
burst period. This p roduction of CNO nuclei w as pre- 
viously discussed bv lHanawa fc Fuiimotol l)1984f K The 
mass fractions at maximum of 14 and 15 are 3.5 % and 
4.3 % respectively. Most of these abundances are large 
because the second combined hydrogen-helium runaway 
is already beginning, but even 1900 s earlier, at 52,668 s, 
the mass fractions were already 0.56 % and 0.99 %, about 
what the CNO processing would give for accreted mat- 
ter with solar metallicity - even though this was a "low" 
metallicity study. 

Because of this, the critical mass for all bursts except 
the first one is smaller, giving a shorter interval between 
bursts. This implies bursts with shorter durations and 
less energy. Less extreme values of density and tempera- 
ture are also achieved and the products of the rp-process 
are not so heavy. Since the second runaway commences 
in the ashes of the first and quickly becomes convective, 
some of the ashes of the first burst, as well as 14 and 15 
catalysts are mixed out into the second accreted layer. 
Since decays have gone on during the inter-pulse period, 
this would be an opportunity for the second burst to 
make even heavier nuclei than the first one. However, 
the lower temperature and density are more critical, and 
the final ashes of the second (and subsequent) pulses are 
actually considerably lighter (Fig. I13|) . 

Because the second burst is more typical for the as- 
sumed accretion rate and metallicity, it is of some inter- 
est to examine its light curve. Fig. El shows that the 
rise to maximum is again relatively slow (compare with 
Fig. occurring on a radiative diffusion time for the 
accreted layer, about 10 s. This is the same behavior 
seen in burst 1 and happens because the convection so 
apparent in Fig. 1101 dies out before the surface luminos- 
ity rises above the background value given by accretion. 
During this convective stage, the star is far from steady 
state. The luminosity at the base of the convective layer 
is orders of magnitude greater than at the top, enabling 
the convec tion zone to extend throughout the layer (e.g., 
I.Tosslll977l) . After this adiabatic condition has been es- 
tablished, the envelope is able to carry the necessary flux 
by radiative diffusion and the rise time for the luminosity 
is slow. As we shall see this contrasts with the situation 



when a hydrogen depleted helium shell runs away f$ l4.2H 
and the luminosity rises almost instantly. 

We shall discuss details of the nuclear physics affecting 
the rise time and tail of the light curve in § 14.1.41 The full 
light curve of the second burst using standard settings is 
given in Fig. 1111 

4.1.3. Bursts numbers 3 through 14 in Model zM 

Bursts 3 through 14 closely resembled burst 2. The 
composition at the onset of the third burst is given in 
Fig. El The burst once again ignites in the outer ashes 
of the previous one. An interesting development is the 
gradual depletion of 12 C which had significant abundance 
in the ashes of the first burst. This is processed to heavier 
elements by a-captures when the heat wave from sub- 
sequent bursts penetrates into the ashes. Fig. and 
1131 show that these subsequent bursts also produce an 
abundance distribution centered on lighter iron group el- 
ements, notably 64 Ge (which decays to 64 Zn), though a 
tail of significant production still extends to A ~ 100. 

Fig.ll4lshows the convective histories for the first three 
flashes, once again emphasizing the difference between 
the first violent, large mass burst and subsequent, mutu- 
ally similar weaker bursts. 

The cumulative bursting history is displayed in Fig.H"5l 
and shows that steady state is achieved after only one 
burst. The composition after 14 bursts is shown in 
Fig. 1161 and the density-temperature structure then is 
given in Fig. El At that time, the only abundances 
lighter than A = 60 having mass fraction greater than 
0.001 are 12 C (0.00 1 6), 28 Si (0.00 68), 32 S (0.0037), and 
60 Ni (0.0052). This very low carbon abundance is insuf- 
ficient to undergo u nstable ignition in deeper lay ers and 
power a superburst l|Cumming fc Bildstenll2001l i. 

4.1.4. Sensitivity of results to nuclear physics at the 
waiting points 

No single reaction rate governs energy generation 
by the rp-process. Early on, energy is produced by 
helium burning and the break out from the beta- 
limited CNO cycle. After the initial destruction of 
14 and 15 0, the burning follows the reaction sequence 
3a ^ 12 C(p, 7 ) 13 N(p, 7 ) 14 0(a,p) 17 F (p l7 ) 18 Nc(a,p) 
21 Na(p,7) 22 Mg(a,p) 25 Al. ... Later, depending upon 
how far up in mass the flow has gone, hydrogen burns at 
a rate sensitive to the capture cross sections, photodisin- 
tegration rates (hence particle separation energies), and 
weak-decay rates for progressively heavier nuclei. The 
rp-process is very similar to the r-process in both its 
dependence upon the properties of waiting point nuclei 
all along the process path, and in that it involves nuclei 
whose properties are poorly determined. 

We undertook a study of the sensitivity of the burst 
light curve to nuclear uncertainties. This is by no 
means the first such study. For other recent efforts see 
Rauscher et "all 120001: iRembeesI Il999t iThielemann et alJ 
20ful Rather than vary a large number of individ- 
ual rates however, we chose, in this initial survey, to 
vary groups of rates. One group is the collection of all 
electron-capture and positron emission rates for nuclei 
heavier than 56 Ni, i.e., the lifetimes of 57 Cu and all heav- 
ier unstable nuclei. This set affects the flow from the iron 
group to heavier nuclei, ultimately near A ~ 100 for the 
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first burst. We varied the rates up and down by one or- 
der of magnitude compared to the standard values, most 
of which were due to measurements or estimates of the 
ground state lifetime. It is likely that our standard val- 
ues overestimate these lifetimes, so multiplying the rates 
by ten is probably more reasonable than dividing them 
by ten, though we did both (Fig. EH- 

It is to be emphasized that the key weak decay rates, 
e.g., for 60 Zn, 64 Ge, 68 Se, etc., are not themselves uncer- 
tain to an order of magnitude, but are probably known, 
even at these high temperatures and densities, to a fac- 
tor of two (§ 12.4(1 . When we vary these decay rates, we 
are really exploring how efficiently the flow goes through 
critical waiting points by a variety of nuclear reactions 
and equilibrium links, e.g., not just Ge{e + v c ) Ga, 
but also, to some extent, 64 Ge(p,7) 65 As(e+i/ e ) 65 Ge, 
64 Ge(p, 7 ) 65 As(p, 7 ) 66 Se (e+u e ) 66 As, etc. The rate 
of these weak flows is sensitive to the quasi-equilibrium 
abundances of nuclei like 65 As and 66 Se as well as their 
half-lives against positron emission. 

Flows may also stagnate before reaching 56 Ni, so wc 
also carried out a set of runs where all rates for unstable 
nuclei heavier than A — 27 were similarly varied. Fi- 
nally, to target more specifically key reaction rates we 
varied individually the rates for groups of nuclei in the 
vicinity of key waiting points around A = 60, 64, and 
68. The results are shown in Fig. EI] In all cases, we 
follow three bursts in model zM. The first burst is more 
directly comparable with the earlier work of lSchatz et alJ 
l(2001af) . while the third may be regarded as more typical. 

Our first observation is that the nuclear rates can af- 
fect the rise time of the burst as well as the tail of the 
light curve (Fig. I2"U|) . In the hydrogen flashes studied 
here, convection dies out before the light curve has risen 
to a small fraction of its maximum. One might expect 
the rise time to then be a consequence of radiative dif- 
fusion. We find, however, that varying the decay rates 
between A = 27 and 56 also has a direct effect on the 
rise. Some key nuclear physics affecting the flow from 
26 Si to A > 56 during the rise of the burst are the de- 
cay rates of 33 Ar, and 3 ' 38 Ca and proton capture on 30 S 
and 34 Ar, both inhibited by photodisintegration of 31 CI 
and 35 K respectively. The flow to 56 Ni passes through 
34 Ar, 35 ~ 38 K, 36 - 39 Ca, 40 - 42 Sc, 44 - 4 3Ti, 44 > 45 V, 45 > 46 Cr, 
47 - 49 Mn, 48 - 50 Fe, 51 - 55 Co, and 52 ~ 56 Ni. 

It is in the peaks and tails of bursts though that the nu- 
clear physics has its most dramatic consequences. Fig. 1181 
and Table show the results for the first burst from 
Model zM. The tail of the burst involves continued hy- 
drogen burning at a rate sensitive to the positron-decay 
lifetimes and proton separation energies of nuclei along 
the rp-process path. After the rise, the values for rates 
between Si and Ni do not appear to be critical, unless 
they arc very small, but different choices for the rates 
above Ni can have dramatic consequences on the shape, 
duration, and peak brightness of the light curve. Most 
important are the flows around A = 60, 64, and 68. The 
leakage through 60 Zn at a relevant temperature and den- 
sity (1.5xl0 9 K; 3xl0 5 gcm" 3 ) proceeds by 60 Zn(p,7) 
61 Ga( P , 7 ) 62 Ge(p, 7 ) 62 Ga(p, 7 ) 63 Ge(e+i/ c ) 63 Ga( P , 7 ) 
64 Ge, with the (p, 7 ) reactions at 60 Zn, 61 Ga and 62 Ga 
strongly hindered by photodisintegration. Most critical 
then are the half- lives of 62 Ge and 63 Ge and proton sep- 
aration energies of 61 Ga, 62 Ge, and 63 Ge. 
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In the vicinity of 64 Ge the critical flows at T = 
1.0xl0 9 K and p = 2.5x10 gem have a different 
character. The dominant reactions are m Gc(e + is c ) 
64 Ga( P , 7 ) 65 Gc(p, 7 ) 66 As(p, 7 ) 67 Se(p,e+ v e ) 67 As(p, 7 ) 
68 Se with critical decays at 64 Ge and 67 Se. Unlike at 
60 Zn, the critical proton captures in this case arc not 
equilibrated, except for 66 As(p, 7 ) 67 Se, and depend on 
cross section more than separation energies. The nucleus 
65 As, long thought to play a c ritical role in admitting rp- 
process flow to heavier nuclei il Wallace fc Wooslevlll 984^1 
is not critical here so long as its proton separation energy 
is as negative as indicated in Tabled The nucleus 66 Se 
is in quasicquilibrium with 65 As at 1.5xl0 9 K but drops 
out by 1.0xl0 9 K so the small abundance of 65 As is a 
hindrance in its production. 

Near 68 Se, the flows are similar to 64 Ge: 68 Se(p,e + v c ) 
68 As(p, 7 ) 69 Se(p, 7 ) 70 Br(p, 7 ) 71 Kr(e+^ ) 71 Br(p, 7 ) 
72 Kr with critical decays at 68 Se and 71 Kr. 

For calculations using faster weak rates, the burst is 
brighter and decays more quickly, as one may expect. 
The converse is true for slower weak rates. Table |S] shows 
that the effects of compositional inertia persist regardless 
of the choice of rates. It is diminished just a little for 
faster weak rates because more hydrogen burns away in 
the outer layers. Conversely it is more dramatic in the 
unlikely case that the rates are much slower. 

The effective lifetime of the waiting point nuclei Ge, 
68 Se, 72 Kr also depends on the rates of two-p roton- 
capture reactions discussed bv lSchatz et alJ 1^.998). The 
reaction formalism applied here treats these processes 
as two sequential proton capture reactions. The re- 
action rates depend very sensitive ly on the associate d 
masses as already pointed out by iBrown et alJ l(2002fl . 
While the Audi-Wapstra masses used here result only in 
a weak flow through the two-proton capture link, the ex- 
perimental uncertainties allow the possibility of a much 
higher reaction flow if 65 As, 69 Br, and 73 Rb are less un- 
bound. The recent di scovery of a + shape isomer in Kr 
l|Bouchez et al.ll2003^ opens the possibility of a substan- 
tial additional reaction flow through such isomeric states 
in the N = Z waiting point nuclei. More experimental 
and theoretical work is needed to investigate that possi- 
bility. 

The termin ation of nuclear flows i n zM are very similar 
to Figure 1 of lSchatz et alJ <|2001a|) including the closed 
loop caused by 107 Te( 7 ,a) 103 Sn. 

4.2. Model Zm 

As an example of a qualitatively different sort of burst, 
we next consider Model Zm. This model had an ac- 
cretion rate five times lower, 3.5xl0 -10 M Q yr _1 , than 
Model zM, and a mctallicity ten times higher (nominally 
"solar" ) . This combination of longer time and a higher 
abundance of CNO catalyst leads to hydrogen depletion 
at the base of the accreted layer long before the first 
burst ignites. Indeed a substantial fraction of the helium 
also burns. At the time of the first burst, 273,780.9 s 
after the onset of accretion, the helium and carbon mass 
fractions at the base of the accreted layer are 30 % and 
67% respectively (Fig. 1221) . 

The runaway commences in the helium shell, near its 
middle for the first burst (Fig. 123(1 , and somewhat higher 
up for the later bursts. Vigorous convection moves out- 
wards colliding after about 0.4 s with the hydrogen shell 
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(—0.32s in Fig. I23C). The mass of the entire accreted 
layer at this point is 6.14xl0 21 g, of which the hydrogen 
layer {X > 0.01) is 7.6xl0 20 g. At this time the high- 
est temperature, 5.0xl0 8 K, is at the base of the helium 
convection zone, where the density is 1.2xl0 6 gcm~ 3 . 
This collision leads to heating if the hydrogen layer and 
an explosion by 12 C(p,7) as (hot) carbon is mixed out 
into hydrogen. Time steps as small as nanoseconds are 
required to follow this interaction. The maximum in en- 
ergy generation shifts from deep in the helium shell to 
the base of the hydrogen shell in which a second connec- 
tive region now grows. Fig. [23 shows this interaction. 
As time passes, the base of the helium zone grows hotter 
at the same time as the hydrogen convective shell digs 
deeper into the helium and carbon. 

The progression of the hydrogen convective shell as 
it follows the helium shell inwards involves some inter- 
esting physics. As the temperature, energy generation 
fFig. I23F ). and luminosity of the helium layer increases, 
it becomes buoyant with respect to the overlaying hy- 
drogen layer and mixing sets in. When the two shells 
connect, briefly mixing 12 C and '"H, rapid energy release 
occurs locally by 12 C(p,7) and a subsequent rp-process. 
This raises the entropy over that of the layers below 
and temporarily shuts off the convective mixing between 
them. But because part of the initial mixing dredges 
down hydrogen, some of this large energy generation oc- 
curs below the original interface, just resolved by one 
or two zones in the present model. With time, burning 
raises the entropy of these zones so that they eventually 
link up with the hydrogen convective shell. The rising 
helium energy generation also keeps the helium convec- 
tion zone in close proximity the the hydrogen shell. In 
this way, the interface moves downward, piece by piece. 
While the specifics of this merger of the two convec- 
tive shells may be model-dependent, it seems inevitable 
that such a merger happens on a very short time scale 
once the helium convection zone encounters the hydro- 
gen layer. The interface becomes, at least episodically, 
Rayleigh- Taylor unstable, resulting in burning, mixing 
"mushrooms". In contrast, the single interface propa- 
gating downward observed in the present calculation is 
due to the limitations of spherical symmetry imposed by 
the one-dimensional Kepler code. A multi-dimensional 
study would be both interesting and, given the short time 
scale, feasible. 

As time passes, the base of the helium zone grows hot- 
ter at the same time as the hydrogen convective shell 
digs deeper into the helium and carbon. The base of the 
hydrogen convective layer also grows hotter and denser 
as the shell becomes more massive. At the same time 
the concentration of heavy elements in the hydrogen 
rises. After 2500 cycles, 70 ms later, the surface lumi- 
nosity climbs to 10 37 ergs _1 . At this point the den- 
sity and temperature at the base of the helium convec- 
tion region, 2.5 x 10 21 g above the bottom of the accreted 
layer, are 3.4xl0 5 gcm -3 and 1.7xl0 9 K. The temper- 
ature and density at the base of the hydrogen convec- 
tive layer, 3.7xl0 21 g above the the iron substrate, are 
2.5xl0 5 gcm -3 and 1.5xl0 9 K. Both shells are generat- 
ing comparable energy. 

One millisecond (600 cycles) later, the helium convec- 
tive shell has been entirely eaten away by the inwardly 
growing hydrogen convective shell. When the the lumi- 



nosity is 10 ergs , the hydrogen shell has moved in 
to 2.5 xlO 21 g, where its temperature and density, at the 
base, are 1.8 xlO 9 K and 2.0 xlO 5 gem -3 . Owing to time- 
dependent convection, there arc gradients in the abun- 
dances in this convective shell, but at the base where 
energy generation is a maximum, the mass fractions of 
1 H, 4 He, 26 Si, and 30 S are 0.036, 0.38, 0.12, and 0.42 
respectively. The inner 2.5xl0 21 g does not participate 
in a major way in the burst. At this point it is still 
mostly helium and carbon. Helium and carbon burn, in 
radiative equilibrium, by an inward moving flame to 28 Si 
fFig. I23|) . Since the fraction of heavier elements ( 28 Si) 
produced by the flame decreases as it proceeds inward, 
the thermohaline instability sets in immediately behind 
the flame. Though this instability is too slow to affect 
the flame itself, it mixes these layers and the even heavier 
ashes of the rp-process above, resulting in a chemically 
homogeneous layer of ash from each burst. The flame 
consumes essentially all of the helium, but ~ 10% of 
carbon remains that continues burning on a longer time 
scale, down to a few percent, as the ashes cool (Fig. I23L S 
and C). 

The heat wave from the convective helium runaway 
also initiates some burning in the non-convective helium 
and carbon beneath before the arrival of the helium burn- 
ing fame. Remaining 14 ' 15 N (from decayed 14 ' 15 0) is first 
converted by (a, 7) reactions to 18 > 19 F which arc con- 
verted by subsequent e + -decay and (0,7) and (a,p) re- 
actions to 22 Ne. This, in turn, is burned into 25 ' 26 Mg 
by another (01,7) or (a,n) reaction. This latter reaction, 
in particular, is responsible for the small "arc" of in- 
creased energy generation visible in Fig. l23D just before 
the helium flame moves in. This arc starts at —0.01 s 
and is caught up by the helium flame at ~ 0.015 s and 
m = 1.75xl0 21 g. 

Seven tenths of a second after the hydrogen and helium 
convective shells first collided (Fig. [531 the temperature 
at the base of the entire accreted layer is 1.75xl0 9 K 
and the density is 7.6xl0 5 g cm~ 3 . Convection is sub- 
siding and most of the hydrogen is gone. The luminos- 
ity is 2.2xl0 38 ergs -1 . Integrated over the entire ac- 
creted layer the dominant abundances are 4 He (0.12), 
12 C (0.055), 16 (0.055), 28 Si (0.14), 34 Ar (0.077), 38 Ca 
(0.15), 39 Ca (0.14), and 40 Ca (0.05). The remainder of 
the burst will be powered by helium burning and, in the 
tail of the burst, the Kclvin-Hclmholtz contraction of the 
cooling ashes. 

After the main burst was over, for example 78 s after its 
onset, the composition consisted chiefly of helium burn- 
ing ashes with appreciable helium itself remaining near 
the outside. Principal mass fractions integrated over the 
accreted layer were 4 He (0.038), 12 C (0.037), 16 (0.032), 
40 Ca (0.27), 28 Si (0.18), and traces of other abundances 
extending up to the iron group (Fig. [23}. 

Subsequent bursts had slightly increasing critical 
masses, larger values of a, but observable properties 
closely resembling those of burst 1. The cause here is 
once again compositional inertia, but of a different sort. 
In the other models, the runaway ignited, to varying 
degrees, in the ashes of the previous burst where the 
composition played a major role in cither catalyzing or 
directly powering the initial nuclear energy generation. 
Here, where ignition always occurs in the freshly ac- 
creted layer, the effect of the accreted ashes upon the 
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opacity of the substrate is more important. Because of 
the lower electron scattering opacity (oc (Z 2 /A)) and 
higher heat capacity (ionic contribution oc 1/ (A)) of the 
ashes of the helium burning layer(s) compared with the 
(assumed) iron substrate at the beginning of accretion, 
the layers below the helium runaway are cooled more ef- 
ficiently in subsequent bursts. Additional helium must 
then burn before the runaway temperature is reached. A 
larger critical mass accumulates to compensate for the 
lower helium mass fraction. At the bottom of the freshly 
accreted layer just prior to the onset of bursts 1-7 the he- 
lium mass fraction was 0.298, 0.266, 0.226, 0.187, 0.150, 
0.137, and 0.105. The remainder was mostly 12 C, some 
(< l%cach) 14 C, 14 > 15 N, 16 < 18 0, 18 < 19 F, 22 Ne and species 
less abundant than lxl0~ 4 . 

This raises the interesting possibility that continued 
burst activity in Model Zm might eventually lead to a 
stable helium burning shell and the accumulation of a 
carbon layer that could power a superburst. Unfortu- 
nately, this was an "expensive" model to follow, requiring 
about 20, 000 cycles per burst and several hundred zones. 
After burst 4 we followed the evolution using a "dczon- 
ing" algorithm that joined neighboring zones when even 
their combined thickness was smaller than 1 cm (mostly 
in the ashes) as they get increasingly compressed by the 
weight of new zones and ashes layers added atop. The 
continued evolution of this model is planned. 

The light curve for the second burst Model Zm is shown 
in Fig. [2U and Fig. Its characteristics, including the 
rapid rise time, were typical for all seven bursts followed. 
The diffusion time of our outer zone (roughly 150 cm at 
10 4 gcm~ 3 ) is approximately 1ms and this sets a lower 
bound on our resolution of the rise of the light curve. 
It seems likely that the rapid rise times seen here will 
characterize all bursts in which a pure helium flash occurs 
beneath a hydrogen layer of appreciable mass. 

Note that in burst 6 the convection reached and stayed 
at the surface, making a more luminous burst (Tabled). 

4.3. Models zm and ZM 

These were similar to Model zM, each being com- 
bined hydrogen-helium runaways. The effects of com- 
positional inertia, though still present were diminished 
relative to Model zM because of the higher initial metal- 
licity (Model ZM; Tabic and longer accretion interval 
(Model zm; Table EJ). The higher metallicity makes the 
critical mass of the freshly accreted layer less sensitive to 
the (nearly solar) metallicity created in the outer ashes 
of the previous burst by primary reactions. A longer- 
burst interval allows the more complete combustion of 
the ashes of the previous burst during the inter-pulse pe- 
riod. 

In fact, the global properties of Models ZM and zM are 
similar in terms of recurrence interval, burst energy, and 
duration emphasizing the greater dependence of outcome 
on accretion rate than composition. After the first flash, 
burst properties are very similar (Fig. I27f) . There are 
differences though. Model ZM has burst intervals about 
10 %-20 % less than zM, but total burst energies (-Eo.oi in 
the tables) 25% smaller. The ratio -Eo.oi /Delay is also 
significantly less. Both these effects reflect the greater 
amount of burning in between pulses for the model with 
the higher metallicity. 

The recurrence times and a values are more regular 



for model ZM than for model zM. This is consistent with 
the fact that the heating in model ZM is determined by 
the accreted metallicity, whereas in model zM, the heat- 
ing is determined by metallicity produced by the previ- 
ous burst. Taam et al. (1993) also concluded that lower 
metallicity led to less regular burst properties. 

The burst energies for Model zm are higher by about 
a factor of three than zM reflecting the larger critical 
mass required for the lower accretion rate. The fraction 
of material that burns between bursts is also greater so 
that the combination Eq.qi/(M x Delay) is smaller. 

Bursts can also be ignited by compositional inertia due 
to helium ignition in the ashes layer. In Model zm the 
energy released by the triple-a reaction in the hydrogen- 
free helium layer of the ashes generated enough energy 
to initiate convective mixing between the ashes layer and 
the overlaying hydrogen layer at runaway. 

4.4. Corrections for General Relativity 

The calculations in thi s paper were carrie d out assum- 
ing Newtonian gravity. iTaam et alJ l(1993t) discuss the 
correctio ns that can be expected for general r elativity 
(see alsolAvasli fc ,TossllT982l iLewin et al.lll993l SlPof 
iCumming et all 120021) For the 1.4 Mq, 10 km neutron 
star we assume here, the gravitational rcdshift correction 
is 1 + z = 1.31. The recurrence time scales and burst du- 
rations in all the tables and figures should be increased 
by (1 + z), the energies decreased by 1 + z, and the lumi- 
nosities decreased by (1 + z) 2 = 1.72. In addition, note 
that the mass co-ordinate used in this paper refers to the 
baryonic mass (number of baryons multiplied by proton 
mass) rather than rest mass. Also the accretion rates we 
give are in baryonic masses in the frame of the simulated 
layer and have to be decreased by 1 + z in the observer 
frame. 

4.5. Mixing by Thermohaline Convection 

Bursts typically leave behind a radial composition gra- 
dient because burning can reach different compositions 
in different zones. For example, a more extensive rp- 
process may occur in outer layers than further in. This 
sometimes happens because of the different density- and 
temperature-dependencies of the 3a reaction, (a,p) reac- 
tions along the rp-process path, and the rp-process itself. 
If "heavier" material - according to its composition - is 
located above "lighter" material, but the compositional 
inversion is less destabilizing than the temperature strat- 
ification, thermohaline convection occurs ( "salt finger 
instability"). Otherwise, if the thermal stratification is 
too week, convection according to the Lcdoux criterion 
sets in. In an ideal gas the compositional buoyancy is de- 
termined by the mean molecular weight gradient ("/«"), 
but in degenerate regions, the mean molecular weight 
per electron ("/i G ") is more important. As a result, when 
the matter becomes more degenerate, thermohaline con- 
vection can switch on or off. Additionally, weak decays 
change /i and /i and their gradients over tim e. We model 
this according to iKinoenhahn et al.l l)1980[) . but using a 
generalization for degenerate equation of state. 

Thermohaline mixing mostly occurs within the ashes 
of each burst, slowly homogenizing its chemical composi- 
tion, frequently only after several subsequent bursts have 
occurred. Eventually, after many bursts, the mixing of 
neighboring layers of ash is also observed. 
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For example, thcrmohaline convection is indicated in 
Fig. 123k for Model Zm. Mixing is first facilitated by the 
behavior of the flame in the helium layer which makes 
less silicon as it moves in. This leads to thermohaline 
convection immediately behind the flame. On the other 
hand, in the hydrogen-rich layer on top, the rp-process 
produces heavier ashes than pure helium burning, leading 
to mixing later on. This later mixing leads to the spikes 
of energy generation after 10 s and between 2x 10 21 g and 
3x 10 21 g in Fig. 116b and are due to mixing events at the 
upper edge of this thermohaline convective region. 

In all other figures showing the convective structure 
we omitted thermohaline convection in the plot in or- 
der to retain the visibility of burning regions, but it was 
included everywhere. 

5. COMPARISON WITH ANALYTIC CALCULATIONS 

We now compare our results with semi-ana l ytic ig - 
nition models, following iCumming fc Bildstenl l)2000|) . 
These models apply a one-zone ignition criterion to sim- 
ple models of the accumulating layer. This approach 
allows a survey of parameter space to be made while still 
giving a good estimate of the ignition depth, and has re- 
cently been applied to observations of th e regular Type 
I bursters 4U 182 0-30 (iCummind 12003]) and GS 1826- 
24 (jGallowav et alJl2003|) . A comparison with the time- 
dependent simulations presented in this paper is of value 
for two reasons. First, it provides a cross-check. Sec- 
ond, it highlights cases in which additional physics not 
included in the semi-analytic models (such as thermal 
and compositional inertia) is needed. 

The models are cons tructed as follows (see 
ICumming fc Bildstenl 120001 for a detailed descrip- 
tion). We integrate the steady-state entropy and heat 
flux equations down through the layer, following the 
composition change as hydrogen burns to helium via 
the hot CNO cycle. Both hot CNO burning and a 
heat flux from below heat the layer. Helium burning 
reactions during accumulation are not included. We 
adjust the thickness of the layer until the "one-zone" 
ignition criterion de^/dT = de rn n\ldT (IBildstenl 119981: 
iFuiimoto et~aTI 119811 : LFushiki fc Larnbl Il987|) islnet at 
the base. Here, e^ a is the triple alpha energy production 
rate (augmented by a factor of 1.9 to account for proton 
captures on carbon), and £ C ooi is a local approximation 
of the radiative cooling. 

The thickness of the layer at ignition sets the burst 
recurrence time. The burst energy is estimated by as- 
suming complete burning of the layer, and a nuclear en- 
ergy release Q nuc = 1.6+4.0 (X) MeV per nucleon, where 
{X} is the mass- weighted mean X in the layer at ignition. 
This expression for Qnuc includes 35 % energy loss in n eu- 
trinos during hydrogen burning flFuiimoto et alj|l987j) , a 
little more than we find in our simulations (see § 12. 5|) . 

Table 1101 compares burst 1 and 2 from each run with 
the analytic predictions. The agreement is very good in 
general. By far the largest discrepancy is for Model zM. 
While the 1st burst agrees well with the analytic predic- 
tions, the recurrence times and burst energies for subse- 
quent bursts are much less than expected. This demon- 
strates very clearly the role of compositional inertia in 
decreasing the recurrence time. As discussed in §0J triple 
alpha reactions in the low density layers during the cool- 
ing tail of the burst give rise to a layer with solar abun- 



dance of CNO (or larger) . Residual hydrogen burning in 
this layer during accumulation gives an additional heat 
source, reducing the mass needed for ignition. 

Compositional inertia is most important for Model zM 
because of the low metallicity. The luminosity 
from deeper layers, Lb, in terms of the energy per 
accreted nucleon, Qb, is Lb ~ 10 34 ergs _1 x 
(Qb/0.1MeV) (M/0.1M Edd ). The hot CNO luminos- 
ity is roughly L CNO ~ £cnoAM h , where £ C no = 
S.SxlO^ergs" 1 x (Z/0M) is the hot CNO energy 
production rate ijFowlerll 196(4 1 Wallace fc Wooslevl ll98lI) 
(Z is the CNO mass fraction), and AMh is the mass 
of the layer containing hydrogen. This gives Lcno ~ 
lO^ergs- 1 x (Z/0.02) (AM H /10 21 g). For the solar 
metallicity models, the hot CNO flux dominates the heat- 
ing of the layer as it accumulates (Lcno ^ -^b), whereas 
for low metallicity, the heat flux from below plays an 
important role in heating the layer (Lcno ~ -^b)- This 
makes the low metallicity models most sensitive to any 
extra heating, e.g., compositional inertia effects. Fig- 
urc r27l shows very clearly the difference between the lumi- 
nosity prior to the 1st burst (~ 10 34 ergs _1 from deeper 
layers) and between subsequent bursts (~ 10 35 ergs _1 
from hot CNO burning of residual hydrogen beneath the 
freshly accreted matter). 

Also apparent in Table 10 is that for Model Zm, the 
Qnuc values are less than expected from burning pure 
helium to iron group (which gives ~ 1.6 MeV per nu- 
cleon). This is due to the fact that much of the he- 
lium burns to carbon prior to ignition. Our analytic 
models do not include helium burning reactions, but 
still give a good estimate of the recurrence time in this 
case. However, they are unable to address the question 
of how much carbon is burned before the runaway oc- 
curs, and whet her a stable helium bur ning layer might 
form. Recently, iNaravan fc Hevll i|2003t) calculated igni- 
tion conditions, including helium burning before ignition 
by approximating the composition profile as that due to 
steady helium burning. They also found that some he- 
lium burned to carbon before instability occurred, and 
referred to these bursts as " delayed bursts" . 

6. CONCLUSIONS 

By coupling nuclear reaction networks of unprece- 
dented size directly to the calculation of zoned stellar 
models of X-ray bursts we have calculated the most re- 
alistic models to date of Type I X-ray bursts on neutron 
stars. Two values of accretion rate and metallicity were 
explored and, in all cases, at least 4 and up to 15 repe- 
titions were followed, assuring that the properties of the 
bursts had reached steady state. The effects of vary- 
ing the nuclear physics were also examined, in particular 
key lifetimes at waiting points along the path of the rp- 
process. 

For the conditions studied, we find, in agreement with 
iSchatz et al.l l)2001aj) . that when the full reaction net- 
work is included, no hydrogen enters the ocean and crust 
of the neutron st ar. Deep hydrogen burning (e.g. as in 
iTaam et al.lll996D . that has also be en discussed as a pos- 
sible energy source for superbursts (|Kuulkers et alJ2002j) 
does not happen and may have been an artifact of the 
limited nuclear reaction networks used in previous stud- 
ies. 

The first burst in each sequence had different prop- 
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erties to subsequent bursts. This is because accretion 
initially proceeds on an inert substrate. However, hy- 
drogen, helium, and CNO nuclei are not depleted in the 
outer layers of each burst and burning continues between 
bursts. This leads to both "the rmal inertia" and "com- 
positional inertia" <|Taamlll98fl . Subsequent runaways 
(Table ignite in the ashes of previous ones, require 
less critical mass and are therefore less energetic. The 
critical mass becomes nearly a constant independent of 
the composition of the accreted matter (at least for mat- 
ter with sub-solar initial composition — see Tableland 
Tables through EJ) . This complicates attempts to in- 
fer the metallic ity of the accreted m aterial from burst 
properties (e.g.. lGallowav et "al"1l2003ft . 

These effects have important implications for the com- 
position of the ashes, which eventually become incor- 
porated into the neutron star crust. In Model zM, 
w hich has the s ame co nditions as the one-zone model 
of iSchatz et al.1 i|2001afl . only the first, energetic, burst 
ends in a SnSbTe cycle. Subsequent bursts ignite sooner, 
reach lower temperatures, and do not produce isotopes 
much heavier than mass 64 (Table 0). The matter that 
accretes into the neutron crust has mean mass in the 60's 
rather than 100. However, it may we ll be that similarly 
violent bursts to ISchatz et alJ l)2001a|) will still be found 
in steady state for different accretion rates. 

As has been pointed out frequently, X-ray bursts are 
marvelous laboratories for the study of nuclear astro- 
physics, especially of nuclei near the proton-drip line 
up to mass 110. Our studies confirm (e.g., § I4.1.4JI the 
sensitivity of the light curves of bursts powered by hy- 
drogen burning to nuclear flows abov e the iron group 
lIKoike et ahlfmM ISchatz et a,1Jl2001aH . We additionally 
find that the rise times of such bursts are also sensitive to 
nuclear decays below the iron group fFig. H8|) . Separating 
these effects out from those due to thermal diffusion and, 
possibly, the spreading of the burning on the neutron 
star will be difficult and will rely on accurate data from 
the nuclear laboratory. However, the light curves for 
Model zM, for example, compare favorably with obser - 
vatio ns of bursts from GS 1826-24 l|Gallowav et all2003fl . 
which Bil dstenl l)2000|) proposed were powered by the rp- 
process. Particularly noteworthy are the long (~ 10 s) 
observed rise times, which compare well with Fig. [201 for 
example. There is much to be learned from a detailed 
comparison of observations and theory for this source. 

In Model Zm (Table d) , the runaway was initiated by 
helium burning beneath a stable hydrogen shell. In con- 
trast to the hydrogen- helium flashes, in which convec- 
tion had ceased by the time the surface luminosity be- 
gan to increase, convection in these bursts continued for 
w 5 seconds after ignition. Also, when the convection 
zone driven by unstable helium burning first broached 
the overlying hydrogen layer, a virtual explosion occurred 
initiated by 12 C(p,7) reactions. The resulting rise time 
was very short, approximately 10 ms. The burst dura- 
tion was shorter than the others, as expected for bursts 
where most of the fuel is helium, and there was a pe- 
riod of super-Eddington luminosity (which would drive 
radius expansion, although this is not followed in our 
models). All four bursts in Model Zm had these same 
characteristics — very rapid rise time, brief duration, and 
super-Eddington luminosities. An interesting question is 
whether these characteristics are shared by all helium 



initiated bursts capped by a layer of accreted hydrogen. 

Our results have implications for carbon-powered su- 
perbursts. For all of the hydrogen-helium flash models 
(zM, zm, ZM), we find very little carbon remains after 
each flash ( < 1% by mass) in a greement with one-zone 
calculations l|Schatz et al.ll20 03). However, in addition, 
we find that carbon is depleted further by alpha-captures 
as the ashes are heated by subsequent bursts. This is be- 
cause some helium remains as well as carbon, and alpha- 
capture is efficient at converting carbon to heavier nu- 
clei such as 28 Si. This implies that unstable hydrogen 
and helium burning at these compositions and accretion 
rates does not lead to accumul ation of sufficient carbon 
to power an unstable r unaway llWooslev fc Taanlll976D 
leading to a superburst l)Cummin^7Tl3iki^to]|l2001[l . In 
addition, the le ss extensive rp-proces s found in Model zM 
compared with ISchatz et alJ ll2001a|) means that photo- 
disintegration of heavy elements jCumming et al.ll2003fl 
will play a less important role in superburst energet- 
ics. Carbon production at these accretion rates may rely 
on some fraction o f stable burning of the accreted fuel 
ISchatz et alJEool . 

In the helium flash model, Model Zm, we find that 
the accumulating helium very nearly burns stably, with 
a substantial amount of helium burning to carbon prior 
to the runaway. This carbon is burned up after helium 
ignition, leaving only ~ 3 % by mass in the deep ashes 
after several bursts. However, the carbon mass fraction 
at the base before the flash was found to increase as 
subsequent bursts were followed, reaching 90% after the 
seventh burst. As we continue to evolve this model, it 
may be that helium burning will ultimately stabilize, 
allowing production of enough carbon to power a su- 
perburst. However, we note that the accretion rate in 
Model Zm is much less than i nferred for the superb urst 
sources (ps 0.1-0.3 Eddington: Ikuulkers et al]l2002ft . 
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Table 1. Proton separation energies of isotones near the long-lived 
waiting point nuclei 64 Ge, 68 Se, 72 Kr, AND 76 Sr. 



Nucleus 5 p (MeV) a Uncertainty(keV) 



65 As 


-0.43 


290 


66 Se 


2.43 


180 


69 Br 


-0.73 


320 


70 Kr 


2.14 


190 


73 Rb 


-0.55 


320 


74 Sr 


1.69 


210 


77y 


-0.23 


unknown 


78 Zr 


1.28 


unknown 



a Taken from the compilation of Brown^^aJJ 120021) . except for the proton separation energies of 77 Y and 78 Zr, which were taken from 
the unpublished calculations of [Brown ( 2 002Th 



Table 2. Summary of Model Properties 



Model 


Z 


Acc Rate 


# bursts 




(Zq) 


(lO" 10 M yr- 1 ) 




zm 


0.05 


3.5 


4 


zM 


0.05 


17.5 


15 


Zm 


1 


3.5 


7 


ZM 


1 


17.5 


12 



Table 3. Composition 81 After Burst 



Model 


Burst num. 


A Z 


X 


A Z 


X 


A Z 


X 


zM 


1 


106 Sn 


0.18 


104 In 


0.09 


106 In 


0.06 


zM 


3 


64 Zn 


0.16 


68 Se 


0.09 


32 g 


0.06 


Zm 


1 


28 Si 


0.27 


40 Ca 


0.23 


24 Mg 


0.11 


Zm 


3 


28 Si 


0.27 


40 Ca 


0.18 


24 Mg 


0.15 



a Three most abundant isotopes at the bottom of the most recently accreted critical mass shortly after the burst indicated is over. 
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Table 4. Flash Properties for Model zM. 



Pulse 


Delay a 


-^pcak 


T .oi b 


Eo.oi 


T l/2 


-El/2 


T ■ C 

' rise 




# 


00 


(lO^ergs- 1 ) 


(s) 


(10 39 erg) 


(s) 


(10 39 erg) 


(s) 




1 


11.59 


1.78 


527 


20.92 


43 


6.10 


2.64 


62 


2 


3.58 


0.89 


382 


6.50 


31 


2.02 


2.24 


61 


3 


3.14 


1.14 


326 


6.37 


25 


2.08 


2.40 


55 


■1 


2.84 


0.97 


308 


5.36 


25 


1.77 


2.49 


59 


5 


2.96 


1.18 


315 


6.23 


20 


1.80 


2.37 


53 


6 


3.21 


0.98 


342 


6.00 


26 


1.85 


2.26 


59 


7 


3.08 


1.02 


339 


6.01 


26 


1.91 


2.50 


57 


8 


2.90 


1.09 


307 


6.03 


23 


1.86 


2.38 


53 


9 


3.46 


0.97 


359 


6.41 


29 


2.02 


2.26 


60 


10 


3.10 


1.13 


330 


6.41 


24 


1.98 


2.38 


54 


11 


3.17 


1.03 


331 


6.11 


24 


1.82 


2.26 


58 


12 


3.17 


1.07 


325 


6.24 


25 


1.95 


2.28 


56 


13 


3.13 


1.02 


339 


6.03 


26 


1.90 


2.37 


58 


14 


3.00 


1.13 


308 


6.22 


24 


1.96 


2.26 


53 


15 


3.52 


1.01 


366 


6.73 


28 


2.02 


2.16 


58 



a "Delay" indicates the time from the last burst or the time since start of accretion for the first burst. L pea ^ gives the maximum luminosity. 
b Duration and total energy for the time during which the burst exceeds 1 % of its peak luminosity are denoted as ro.oi and Eo.Oi- The 
time and energy for emission above 50% maximum are tj/2 an d i?i/2 
c Timc for the burst to rise from 10 % to 50 % of its peak luminosity. 

d Persistent luminosity compared to integrated burst energy (here we integrate over the time the burst has 1 % of it peak luminosity). 



Table 5. Model zM with weak rates (A > 56) ten times faster. 



Pulse Delay 


^peak 


T0.01 


-Eo.oi 


Tl/2 


£l/2 


Trise 


Q: 


# 


00 


(10 38 ergs 


- 1 ) (s) 


(10 39 erf 


E) (b) 


(10 39 erg) 


(s) 




1 


11.60 


2.30 


332 


20.81 


82 


15.40 


3.42 


62 


2 


3.43 


0.95 


261 


6.25 


53 


3.66 


2.02 


61 


3 


4.07 


1.24 


251 


7.92 


55 


5.14 


2.09 


57 


4 


4.00 


1.19 


251 


7.65 


55 


4.85 


2.09 


58 




Table 6. 


Model zM with weak rates (A > 56) ten times slower. 




Pulse Delay 


^peak 


TO. 01 


Eo.Ol 


r l/2 




''rise 


ct 


# 


00 


(10 38 ergs 


- 1 ) (s) 


(10 39 erf 


;) (s) 


(10 39 erg) 


(■) 




1 


11.59 


1.75 


923 


15.84 


44 


5.94 


2.60 


81 


2 


2.75 


1.10 


579 


10.11 


31 


2.38 


1.44 


30 


3 


2.70 


1.32 


357 


6.14 


20 


2.02 


0.82 


49 


4 


6.82 


2.30 


424 


10.31 


16 


2.88 


2.22 


74 






Table 7. 


Flash Properties for Model Zm. 






Pulse 


Delay 


-^pcak 


TO. 01 


-Bo. 01 


Tl/2 


#1/2 


T ' a 

' rise 


O: 


# 


00 


(10 38 ergs -1 


) (b) 


(10 39 erg) 


(s) 


(10 39 erg) 


(ms) 




1 


76.05 


3.33 


103 


7.37 


28 


6.73 


32.07 b 


231 


2 


80.43 


4.73 


83 


7.64 


13 


6.06 


1.38 


235 


3 


86.14 


5.06 


87 


7.95 


12 


6.19 


0.91 


242 


■1 


96.19 


5.10 


93 


8.07 


12 


6.14 


1.04 


267 


5 


107.00 


4.55 


103 


8.54 


15 


6.86 


1.82 


280 


6 


111.20 


6.33 


27 


7.98 


11 


6.52 


0.61 


312 


7 


130.84 


5.12 


107 


9.20 


14 


6.87 


0.51 


318 



a Note that the rise time in this table is in ms not s as in the other tables. 

b This first burst shows a triple peak. If only the time required to rise to 50 % of the first peak is considered, the corresponding rise time 
is less than 0.7 ms. Indeed, if the rise time to 30% of peak is considered, all bursts have rise times between 0.3 ms and 0.4 ms, except for 
the first burst that rises in 0.7 ms, as mentioned above. 
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Table 8. Flash Properties for Model ZM. 



Pulse 


Delay 


^peak 


TO. 01 


Eo.oi 


T l/2 


El/2 


Trise 


a 


# 


(h) (lO^ergs -1 ) 


(s) (10 39 erg) 


(s) 


(10 39 erg) 


(s) 




1 


3.47 


1.91 


124 


4.93 


15 


2.15 


0.66 


78 


2 


2.66 


1.51 


141 


4.61 


16 


1.88 


0.55 


64 


3 


2.72 


1.58 


143 


4.65 


17 


2.01 


0.56 


65 


1 


2.69 


1.51 


138 


4.59 


18 


2.01 


0.59 


65 


5 


2.65 


1.63 


133 


4.68 


15 


1.94 


0.54 


63 


6 


2.74 


1.54 


135 


4.62 


17 


2.02 


0.58 


66 


7 


2.65 


1.52 


143 


4.60 


16 


1.87 


0.59 


64 


8 


2.70 


1.50 


146 


4.66 


17 


1.96 


0.57 


64 


9 


2.65 


1.57 


144 


4.62 


17 


2.02 


0.57 


64 


10 


2.68 


1.55 


141 


4.58 


16 


1.92 


0.55 


65 


11 


2.68 


1.65 


134 


4.72 


17 


2.09 


0.51 


63 


12 


2.73 


1.66 


133 


4.70 


16 


2.04 


0.51 


65 




Table 9. Flash Properties for Model 


ZM. 






Pulse 


Delay 


^peak 


TO. 01 


fio.oi 


Tl/2 


El/2 


Trise 


a 


# 


(h) (lO^ergs -1 ) 


(s) (10 39 erg) 


00 


(10 39 erg) 


(s) 




1 


60.16 


2.97 


313 


20.92 


37 


8.23 


1.09 


64 


2 


57.27 


2.36 


333 


20.43 


55 


11.14 


0.87 


62 


3 


54.91 


2.61 


301 


19.48 


48 


10.20 


0.93 


63 


4 


53.03 


3.13 


275 


18.88 


36 


8.44 


1.08 


(-"; o 


Table 10. Comparison of burst properties with analytic models 








Trecur 


Energy 


Qnuc 


(X) b 


(Y) b 










00 


(10 39 erg) 


















Model zm 














burst 1 


60.2 


20.9 


4.51 












burst 3 


54.9 


19.5 


4.60 












analytic 


78.0 


24.5 


4.09 


0.623 


0.376 












Model Zm 














burst 1 


76.1 


7.37 


1.26 












burst 3 


86.1 


7.64 


1.20 












analytic 


70.3 


9.73 


1.80 


0.051 


0.929 












Model zM 














burst 1 


11.6 


20.9 


4.68 












burst 3 


3.14 


6.37 


5.25 












analytic 


14.0 


24.3 


4.54 


0.735 


0.264 












Model ZM 














burst 1 


3.47 


4.93 


3.68 












burst 3 


2.72 


4.65 


4.43 












analytic 


4.64 


6.72 


3.77 


0.544 


0.436 







a The nuclear energy release in MeV per accreted nucleon. For the analytic models, we set Qnuc = 1.6 + 4 (X) McV. For the time-dependent 
models, Q nuc = E . 01 /(t iccui M) = 2.59 MeV/nucleon X (£ . i/10 39 erg) (r reC ur/l h^^M/l^ X lO" 9 M yr" 1 )- 1 . 
b The mass-weighted hydrogen or helium fraction immediately before the burst, e.g., (X) = i J X(M)dM . 
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Fig. 1. — Color-coded nuclear mass excesses taken from different data sources in the region of interest for the rp-proccss. A dar k Jine 
shows the location of nuclei with Z = N. The actual size of the networ k used in these studies varied with time and model 12,11 , Red 
triangles indicate nuclei with experimentally determined mass excesses (Audi^i Wapstra 1993). Nuclei with black circles surrounding a 
symbol indicate that an extrapolated or interpolated mass is available from TATichT^W S^^^ ^ For those ci rcles enclosing a re d triangle, 
these estimated values were used. Green squares, circled or not, show where data from TBrowri et alj i2002l) and IBrownl (2002) was used. 
These add a calculated displacement energy to the Audi- Wap stra masses for N > Z to obtain the masses of mirror nuclei with Z > N. 
For a subset of these (black squares inside green squares) from lBrown et "all Ooo3) the errors in the calculated displacement energies are 
~ 100 kcV. The soJid green squares which cover a wider region are from lBrowri 120021) and have a larger error of several hundred keV in the 
calculated displacement energy due to the wider extrapolation and the possible effects of deforma tion in the A = 8 mass region. In the 
region of N = Z from A = 76 — 100, the error in the mass excess is dominated by the errorin_the|Aud^&; Wapstra extrapolations, which 
arc on the order of 0.5 MeV. Finally, bJue circles indicate mass excess data taken from ICToller et al ■I 119951) . These were used wherever 
shown, (circled or uncirclcd). 
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Fig. 2. — Systematics of weak decays for nuclei along the rp-process path near 72 Kr. Nuclei labeled "(Fermi)" in this figure are those 
whose decay rates are dominated by Fermi transitions and expected to be insensitive to the temperature. Nuclei labeled "(Fermi-P)" are 
those whose decay rates decrease quickly with temperature. Estimates of the low-temperature decay rates for these nuclei depend largely 
on estimates of the nuclear partition functions. Nuclei labeled "(2+)" are those even-even nuclei whose low-temperature decay can be 
sensitive to the decay rate of the first excited 2+ state. The number in each box is the estimated ratio of the total (e+ + electron capture) 



decay rate at T = 1.5xl0 9 K, p = 
nuclear flows during the rp-proccss. 



10 g cm 14 to the ground state e+ decay rate. The dashed black line illustrates important weak and 
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Fig. 3. — Composition just prior to the first burst of Model zM (at 41, 640.5 s). Combined hydrogen burning by the beta-limited CNO- 
cycle and helium burning are responsible for a rapidly rising energy production. The total nuclear power here is 3.4xl0 35 ergs -1 and the 
time before the burst, about one minute. 
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Fig. 4. — The first burst in Model zM. The lower panel shows that ignition initially occurs above the composition interface, but grows 
over the next 6 s, both outwards and inwards, to encompass the entire shell. Green hatched regions are convective; red cross hatched 
regions are semi-convective. The observed burst commences shortly after convection reaches the surface and has begun to recede so that 
the luminosity during the rise is chiefly transported by diffusion. Times are offset by 41, 700 s since the beginning of accretion. Any effects 
due to the spreading of burning over the surface of the neutron star are ign ored . Though given here for our standard choice, the rise time 
is sensitive to the nuclear reaction rates empl oyed in the calculation (Fig. |2HJ. In this and all subsequent depictions of the light curve, 
general relativistic effects have been ignored I - !; 14.11 . 
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Fig. 5. — The top figure shows the light curve of burst 1 of Model zM. Following the main pulse, lasting about a minute, there is a long 
tail, lasting perhaps 5 min, powered by the continued burning of hy drogen. The shape of the light curve and its tail are sensitive to nuclear 
decay rates and proton capture along the rp- process path ("Fig. 1181 , The bottom panel shows that appreciable convection only occurs prior 
to the rise of the pulse (see also Fig. |1J . Shades of blue color indicate the nuclear energy generation rate while shades of purple indicate 
energy loss by neutrinos, both on a logarithmic scale. The flow of heat into the neutron star substrate is apparent. 
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Fig. 6. — Composition and nuclear energy generation during the convective phase of burst 1 in Model zM (41,705.4 s). The rp-process is 
responsible for producing energy and intermediate mass nuclei, the most abundant two of which are here 18 Ne and 24 Si. Nuclei as heavy 
as 64 Ge already have appreciable abundance (2xl0 — 4 by mass). The region of near constant abundances is convective. Nuclear power 
is 8.3xl0 37 erg s —1 , though the luminosity transported to the surface is still only 1.3xl0 35 ergs - 1 , invisible compared with that due to 
accretion. The observable burst has not started yet. 
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Fig. 7. — Composition of Model zM and nuclear energy generation in the tail of the light curve of burst number 1 (41, 992 erg s _1 , Fig.l5l, 
Though depleted at the base, considerable hydrogen remains at higher altitudes and burning continues at a high rate. The ashes at the 
bottom of the layer are rich in 106 Sn, the termination of the rp- process. 
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accreted mass / 10 g 

Fig. 9. — Composition and energy generation at the onset of the second burst in Model zM (5,4592 s), 29 s before the flash. The ignition 
of the second burst occurs in the ashes of the previous one, not at the bottom of the newly accreted layer (which is located at 4.7x 10 21 g). 
This occurrence makes the critical mass for second and subsequent bursts smaller and less sensitive to the metallicity of the accreted 
material. Temperature and density at the location of maximum energy generation are 9.4x10 s gem -3 and 2.9x10 s K respectively. A short 
time later, convection mixes the heavy ashes, e.g., 64 Zn, out into the freshly accreted layer. 
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Fig. 10. — Rise time and convection for burst number 2 of Model zM. This is a more typical burst for the model than the one shown 
in Fig. [1] Once again, convection has ceased by the time the burst first becomes visible. Following a brief convective stage lasting about 
2 s, well above the bottom of the freshly accreted layer, the principal burning starts inside the ashes of the previous burst which end at 
4.7xl0 21 g (Fig.|5J. The intensity of the blue scale indicates nuclear burning. Note that appreciable burning occurs in the ashes of burst 
number 1 as the heat wave from burst number 2 propagates through. This behavior is in contrast to th e ne utrino losses which dominated 
in Fig. |2] The rise time is sensitive to the nuclear reaction rates as well as the diffusion time scale (Tig. 120^ 



2G 



Woosley et al. 




2 - 



- 1 1 I 



convection 
0000 semiconv. 



100 



200 



300 



t/s 



54600 



< 



to 
m 



o 



; 



C/) 

to 
o 



Fig. 11. — Light curve for burst number 2. This briefer, less luminous second burst is more typical of all subsequent bursts in Model zM. 
Shades of blue indicate nuclear burning. 
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Fig. 12. — Mass fractions and energy generation at the onset of the third burst in Model zM. This shows that the conditions at which 
the third burst ignites are very similar to those of the second. Compare with Fig. |5] 
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and 104 Pd. In burst 3 they were 64 Zn and 68 Ge. Abundances are evaluated at the base of each accreted layer after decay. 
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"~ time *~ 

Fig. 14. — Kippenhahn diagram for bursts 1, 2, and 3 in Model zM. Green hatching indicates convection; red cross hatching semi- 
convective regions; blue shading shows net energy generation (nuclear minus neutrino losses); and pink shading indicates net energy losses 
(neutrino losses exceeding nuclear energy generation). Each level of blue color indicates an increase by one order of magnitude. The y-axis 
gives the enclosed mass coordinate above the assumed neutron star substrate and the thick black line gives the total mass of accreted 
material. The x-axis gives time increasing from left to right, but different parts of the evolution are plotted on different time scale. Breaks 
on the axis indicate a change of time scale and and below each segment we give the length of that time interval (not the total time). The 
rate of change of the total mass is inversely proportional to the magnification of time in each section and is thus also and indicator of the 
evolutionary time scale in each section. Note repeated waves of nuclear burning in the ashes of previous bursts as heat from the current 
burst propagates inwards. There are also periods of convection in between bursts. 
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Fig. 15. — Fourteen flashes from Model zM. Note the regularity of the last 13. The first is a start up transient. Note also the heat flow 
and burning in the ashes of earlier bursts. 
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accreted mass / 10 g 

Fig. 16. — Composition of the entire accreted shell after the 14 th burst in Model zM. This is the composition that will ultimately be 
merged into the neutron star. Aside from the anomalously heavy ashes associated with the first burst, that composition is mostly ° 4 Zn 
and 68 Ge. Very little unburned carbon remains in these ashes. 
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Fig. 18. — Sensitivity of the light curve of the first pulse in Model zM to five choices of weak rates - standard, standard times 10 above 
A = 56, standard divided by 10 above A = 56 and similar modifications above A = 27. Though only the weak rates were altered, the 
changes reflect uncertainty in all processes, including proton capture inhibited by photodisintcgration, that affect the nuclear flow. The 
shape of the light curve is clearly quite sensitive to the nuclear data set employed. 
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Fig. 19. — Sensitivity of the light curve of the first pulse in Model zMto variations along the waiting points in the vicinity of A = 60, 
64, and 68. The nominal light curve is shown along with the result when all weak rates above A = 59 are multiplied by 10 (sec Fig. 1181 
changing A m i n from 57 to 59 has no effect). Also shown are the results of progressively adding in accelerations to flows in the mass ranges 
A = 60-63, 64-67, 68-71, and 72—79. Details of the flows for these mass ranges are given in the text. A separate calculation, not shown, in 
which only the decay rate of 64 Ge was accelerated by 10 is virtually indistinguishable from the curve e"*~, EC(59 < A < 68) xlO. Factors 
affecting leakage out of this single nucleus thus dominate the flow from the iron group to 68 Se. The "blip" at 180 s is the addition of a new 
surface zone by accretion. 
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Fig. 20. — Sensitivity of the rise of the light curve of the third pulse in Model zMfor three choices of weak rates - standard, standard 
times 10 above A = 27, and standard divided by 10 above A = 28. Flows affecting the rise time are discussed in the text. Time zero is 
denned as when each burst reaches 10 37 ergs — 1 . Because this is the third burst, there have been cumulative effects from the altered rates; 
the critical masses of the burning layers, for example, are not the same (Tables HI- 151. nor are the total burst energies. 
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Fig. 21. — Sensitivity of the light curve of the third pulse in Model zMto variation in the nuclear physics employed (see Fig. 1181 . This 
third burst may be more representative than the first one. 
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Fig. 22. — Composition at the beginning of the first burst in Model Zm. Runaway here occurs in the helium burning shell where a 
lot of carbon (67 % by mass) has already been synthesized. Subsequent convective mixing of the helium and hydrogen layers leads to the 
explosive burning of hydrogen in the presence of a large amount of carbon catalyst. The final products will not be particularly heavy. 
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Fig. 23. — Convective structure of the first burst generated in Model Zm viewed on progressively smaller time scales. Note the adjustment 
of the definition of zero time in each frame. A: On a time scale of minutes the burst is followed by an extended period of slow burning 
(indica ted by blue color) and mixing by thcrmohalinc convection (yellow cross-hatching); B: The light curve itself lasts about 13 s (similar 
to Fig. 1261 which is for the second burst) and the burning layer is at least partly convective about half that time. The runaway and 
convection both begin well above the bottom of the accreted layer. The inner 2.3xl0 21 g of carbon-rich material is never convective, but 
does burn radiatively (note the dark blue band between 2 and 4.5 s). C: At still higher resolution, one sees the growth, both inwards 
and outwards, of the helium burning shell and its collision with the hydrogen layer at —0.32 s at 5.38xl0 21 g. The collision causes a mild 
explosion and rapid growth in the convection zone which bifurcates into hydrogen-rich and helium convective shells. D-F: Subsequent 
frames, at increasing resolution, show the inward propagation of the base of the hydrogen convective shell. Many thousands of stellar 
models were needed to follow this in this one burst. By the end, the hydrogen and helium shells are completely merged. In the radiative 
helium-carbon layer, at about —0.01s on Panel D, a weak front of 22 Ne burning by 22 Ne(cs,n) and 22 Ne(cr,7) begins to move inward. 
This is followed by a self-sustaining helium burning flame ignited at +0.01 s. This flame leaves behind mostly silicon and carbon, but little 
helium. Carbon continues to burn after the flame has passed, mostly between —0.28 s and —0.05 s on Panel C. The fact that convection 
persists (above m > 2.5xl0 21 g), carrying high luminosity almost to the surface for over 0.1s leads to a very rapid rise in the light curve. 
Unlike Model zM, where compositional "inertia" plays a big role, the critical masses, convective structures, compositions, and light curves 
are very similar in subsequent bursts from Model Zm. 
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Fig. 24. — Composition following the first burst in Model Zm. The time is 78 s after the beginning of the burst, well into the tail of the 
event whose brilliant display lasted about 30 s. Some slow burning is still occurring in the helium layer. Interior to 4xl0 21 g the layer is 
well mixed by thermohalinc convection due to a mild inverted gradient in A. Only 4 % carbon remains in these inner ashes, but this will 
survive, since there is no helium there. Subsequent bursts leave similar ashes. 
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Fig. 25. — Rise of the light curve for the second burst generated in Model Zm. The light curve of the first burst is similar. The rise from 
10 37 erg s -1 to X 10 38 ergs takes less than 1 ms (neglecting the propagation time around the neutron star). The rise from lx 10 38 erg s — 1 to 
2xl0 38 ergs -1 takes about 100ms. Above the Eddington luminosity. 2xl0 38 ergs -1 , there will be radius expansion not properly followed 
in the present study. The excess energy will go into driving a wind. 




Fig. 27. — Repeated flashes from Model ZM. During the inter-pulse period the luminosity from accretion (2.0 X 10 37 erg s 1 ) would 
obscure the thermal emission from the cooling ashes that is shown. 



